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SPACE  PROTON  DOSES  AT  POINTS  WITHIN  THE  HUMAN  BODY 


David  L.  Dye 
The  Boeing  Company 

Abstract 

Distribution  patterns  of  absorbed  dose  influence  the 
radiation  response  of  a  mammalian  system,  because  of  differing 
radiosensitivities  of  different  organ  systems.  In  a  man  ex¬ 
posed  to  space  radiations;  e.g. ,  an  astronaut,  body  self¬ 
shielding  produces  nonuniform  dose  distributions  which  depend 
upon  external  shielding  configurations  and  the  Incident  radi¬ 
ation  parameters.  In  this  paper  are  presented  the  doses  at 
twelve  specific  points  in  the  body  of  a  seated  man  exposed  * 

to  isotropic  incident  space  protons,  where  the  man  is  inside 
* 

various  thicknesses  of  external  vehicle  shell  shielding.  The 
body  points,  selected  for  their  radiobiological  interest,  are 
in  (or  on)  sternum,  chest  skin,  femur,  spinal  column,  eye, 
central  gut,  and  a  series  at  various  lateral  depths  on  the 
waist.  The  protons  reaching  these  specific  points  from  all 
directions  traverse  tissue  thicknesses  that  were  determined 
from  scale  drawings  of  a  statistically  standard  man  (75-p«rcea- 
tile).  The  proton  penetration,  secondary  radiation  generation, 
and  total  dose  delivered  to  each  specific  body  point  was 
calculated  using  an  .IBM  (Fortran)  computer  code. 

For  feasibly  thin  aluminum  shields  and  typical  spectra,  the 
secondary  dose  is  a  small  fraction  of  the  total  proton  dose  for 
most  of  the  body  points.  The  doses  at  a  given  point  vary  with  both 
spectral  shape  and  external  shielding,  as  well  as  point  to  point. 
Power  law  spectra,  S  adE,  normalised  to  one  proton/cm^-sec  of  E  >  100 
Mev,  with  2  ^  a  4  5»  sad  the  Freden  and  White  trapped  proton 
spectrum,  were  considered.  Resulting  dose  values  for  n  ■  3  to  4  (a 
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typical  solar  flare-associated  proton  spectrum)  vary  from 
—4  p 

about  10  to  10  rad/proton/cm  >  100  Mev  for  the  various 
body  points  for  a  1  gm/cm2  external  shield,  and  from 
10-7  to  10  rad/proton/cn2  >  100  Mev  for  a  10  gm/cm2 
external  shield. 


Introduction 

Dose  distribution  patterns  in  a  biological  system  have  an  important 
influence  on  the  radiation  effects  in  that  system.  In  the  human  body 
exposed  to  a  typical  space  proton  environment,  there  is  some  self¬ 
shielding  that  reduces  the  body  interior  doses  with  respect  to  the 
skin  doses  and  thus  affords  some  measure  of  radiation  protection.  Pre¬ 
vious  studies  have  been  made  of  the  depth-dose  patterns  resulting  from 
several  simple  proton  spectra  impinging  on  slab  shields  and  penetrating 
into  tissue  slabs.  ’  These  early  results  show  the  dependence  of  depth- 
dose  on  spectral  shape  and  exterior  shielding,  but  the  simplified 
geometry  used  leaves  one  still  wondering  just  what  the  absorbed  doses 
are  inside  a  biologically  interesting  shape  such  as  the  human  body. 

The  earlier  results  were  calculated  by  means  of  an  IBM  7090 
computer  code.2  A  slightly  modified  version  of  this  same  program  was 
used  in  part  of  the  present  calculation.  A  principal  feature  of  the 
present  program  is  that  it  calculates  a  series  of  dose  values  due  to 
penetrating  protons  in  a  set  of  slab  geometries,  then  sums  these  dose 
values,  weighted  as  neoeseary  by  the  fractional  solid  angles  subtended 
by  eaoh  slab  problem,  to  obtain  the  proper  total  dose  at  an  interior 
point.  Tor  the  work  here  reported,  the  seated  human  (75-percentile 
male)  body  was  used,  with  various  spherical  shells  around  it.  Doses 
were  computed  at  a  number  of  body-interior  points:  femur,  sternum. 


1.  '  D.  L.‘  Dye  and  J.  C.  Noyes,  Biological  Shielding  for  Radiation 

Belt  Particles,  J.  Astron.  Sci.  VII,  64  (1959) • 

2.  D.  L.  Dye  and  0.  Butler,  Computer  Calculations  of  Doses  from 
Protons  in  Spaoe,  £.  Astron.  Sci.  IX,  63-71  (19&)* 
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backbone,  eye  lens,  central  gut,  and  at  waist  level. ^  The  results  of 
these  computations  are  tabulated  here  and  show  a  dependence  on  spectral 
and  shielding  parameters  similar  to  the  earlier  results.  However,  these 
values  now  represent  useful  doaea  (per  isotropic  incident  proton/om  of 
E  >  100  Mev)  at  the  specific  body  points.  If  the  incident  isotropic 
flux  and  spectrum  is  known  exactly— an  unlikely  eventuality— then  the 
proton  ionization  dose  at  any  of  these  points  may  be  known  accurately. 

The  dose  due  to  secondaries  generated  in  nuclear  interactions  has 
been  estimated  variously  but  it  is  small  for  thin  shields.  A  com¬ 
puter  program  has  been  developed  to  determine  the  tissue  doses  due  to 
the  secondaries  generated  by  the  incident  primary  protons  in  the  shield¬ 
ing  material.  The  secondary  component  includes  inelastic  scattering- 
produced  gamma  rays,  low-enersy  evaporation  nucleons,  and  high-energy 
knock-on  nucleons.  Part  of  the  tertiary  neutron  and  gamma  effects  were 
included,  by  adding  the  secondary  proton  spectrum  to  the  primary  at 
every  point,  so  the  total  was  used  in  interior  sublayers.  The  results  • 
Justify  the  usual  procedure  of  neglecting  the  secondaries  from  incident 
space  protons  of  typical  energy  spectra. 


The  Computer  Programs  ' 

Primary  Dose  Code 

Spectral  Modification.  In  the  present  version  of  the  primary 
proton  dose  calculation,  an  incident  proton  spectrum  is  degraded  through 
shielding  layers  according  to  formulas  1  to  3  of  Ref.  2.  However,  in¬ 
stead  of  interpolating  between  tabulated  range-energy  points,  as  was 


3.  D.  L.  Dye,  A  Geometrical  Analysis  of  the  Seated  Human  Body  for 

Use  in  Radiation  Dosage  Calculations.  The  Boeing  Company,  D2-9010?, 

T19657. 

4.  R.  I.  Allen,  et  al.,  Shielding  Problems  in  Manned  Space  Vehicles. 
Lockheed  Aircraft  Corp.,  NT-1 40,  (1961). 

5.  R.  K.  Wilson,  et  al.,  A  Study  of  Space  Radiation  Shielding  Problems 
for  Manned  Vehicles.  General  Dynsnios  Corp. ,  FZK-lMt  (1962). 


-4- 


done  in  Ref.  2,  we  have  made  a  least  squares  fit  to  those  tables  to 
obtain  analytic  expressions  for  proton  range  as  a  function  of  energy  and 
energy  as  a  function  of  range.  Thus,  the  proton  range  function  is  de¬ 
fined  empirically  to  bes 

RNGE?(2)  =  exp(Al  +  A2  log  E  +  A3  log2  E),  (1) 

and  the  proton  energy  function  is  (by  quadratic  inversion): 

EOT®  Sia  -  v»B ill,  «> 

where  E  is  proton  energy,  R  is  proton  range,  and  the  A's  are  coefficients 
of  the  fit  which  depend  on  the  shielding  material*  Equation  1  was  chosen 
as  the  form  to  which  to  fit  by  least  squares  the  range-energy  tables  be¬ 
cause  range  is  approximately  a  power  law  function  of  energy;  that  is, 

R  ■  RE0,  (3) 

with  k  and  n  slowly  varying  functions  of  E  over  the  range  of  interest, 

1  £  E  ^  1,000  Mev. 

Use  of  an  analytic  range  function  provides  simply  for  the  'solution 
to  the  spectral  modification  with  penetration  into  shielding  material. 

It  also  allows  one  readily  to  compute  the  output  spectrum  on  a  speci¬ 
fied  energy  grid,  so  that  spectra  are  easily  added  at  an  interior 
point. 

If  a  monodirectional  proton  speotrum  P(E)dE  is  Incident  on  a  shield 
layer  of  thickness  X,  then  the  penetrating  speotrum  P'(E')dE'  may  be 
computed  by  means  of  the  recipe: 

E'  -  ENGYF(RNGEF(E)  -  X),  (4) 

P'(E')dE»  -  P(E)dE(d£'/dE),  (5) 

and  (dR/dE)_ 

(dS'/dS)  a  (6) 


6.  M.  Rich  and  R.  Madey,  Range  Energy  Tables.  UCSL-2301  (1954) 


-5- 


where  the  subscripts  mean  that  the  derivative  is  evaluated  at  E  and  E' , 
respectively. 

Zn  the  present  ZSM  code,  an  output  energy  grid  S' 00  ia  speolfied, 
so  the  E(K)  values  needed  to  obtain  it  are  computed  by  Eq.  4.  Then  the 
incident  specified  spectrum  P(K)  is  rearrayed,  by  interpolation,  onto 
this  input  energy  grid.  Then  the  "energy  spread  factor",  determined 
by  Eq.  6,  is  multiplied  into  the  spectrum  P  to  obtain  the  output  spec¬ 
trum  ?'  arrayed  on  the  specified  output  energy  grid  E'(K),  according 
to  Eq.  5. 

This  spectrum  degrading  process  is  repeated  for  each  layer  of 
shielding  material  until  the  point  of  interest  for  dose  calculation  is 
reached.  The  actual  dose  computation  is  done  as  it  was  earlier  in 
versions  of  this  code;  namely,  by: 

/  1,000  Mev 

P'(E')  DOSF(E')  dE*  (?) 

1 

where  DOSF(E)  is  proportional  to  d^dX  for  tissue. 

Vehicle  Analysis.  It  is  desired  to  compute  a  dose  absorbed  by 
some  material  such  as  tissue,  within  a  complicated  arbitrary  ehielding 
configuration  such  as  a  space  vehicle  or  a  human  body,  due  to  an  arbi¬ 
trary  angular  dietribution  of  incident  proton*.  The  way  this  is  done 
may  be  seen  by  referring  to  Fig.  1.  Briefly,  the  three-dimensional 
material  configuration  ia  analysed  in  a  large  number  of  one-dimensional 
problems,  so  that  the  oode  just  described  can  solve  the  three-dimensional 
problems. 

Looking  outward  from  an  interior  .point  ?  in  all  directions,  on* 
sets  various  solid  angular  regions  AA,  ovsr  each  of  which  th*  shield¬ 
ing  layers  are  approximately  the  same.  If  the  radiation  reaching  P  ia 
direct,  or  primary  and  unsoattered,  then  we  may  approximate  the  dose 
at  P  by  the  sum: 


6 


F;GCS2  1 


TYPICAL  SPACE  VEHICLE  SHIELDING  THICKNESS  ANALYSIS 


(8) 


DOSZ(?)  =  F  [l)OSS(SIp)]  [aNG(IP)]  [^Xp] » 

I?=l 

where:  DOSE(S^p)  is  the  absorbed  dose  computed  behind  a  layer  shield  Sjp» 

represents  the  one-dimensional  shielding  configuration 
of  whatever  layers,  materials,  and  thicknesses  speci¬ 
fied  for  the  solid  angular  region  indexed  by  IP, 

ANG(IP)  is  the  fraction  of  the  flux  coming  from  the  IP  direction 
calculable  from  the  angular  distribution  of  the  incident 
radiation,  and  normalized  to  be  a  constant  unity  value 
for  isotropic  flux, 

F  is  the  incident  flux  assumed  suitably  no:  Aized, 

is  the  fractional  solid  angle  subtended  about  the  point 
of  interest  P  by  the  constant  shielding  region  IP,  and 

NP  is  the  total  number  of  such  regions  needed  to  complete 
the  full  solid  angle,  so  that  the  vehicle  shielding 
is  analyzed  into  sufficiently  small  angular  regions 
to  obtain  the  desired  accuracy. 

In  Fig.  1,  six  regions  were  chosen  to  illustrate  the  concepts,  but 
an  actual  analysis  might  use  more.  In  the  body  point  results  reported 
here,  1$  to  25  regions  were  used  in  most  cases.  It  was  found  that 
little  change  in  the  final  dose  answer  is  obtained  with  larger  numbers 
of  subregions  than  this;  that  is,  the  dose  calculation  converges  with 
relatively  large  Ail  values  for  practical  cases.*  Equation  8  describes 
the  method  of  calculating  a  dose  to  an  interior  point  due  to  incident 
primary,  straight-through,  radiation  of  arbitrary  incident  angular 


*See  the  paper  by  J.  W.  Keller,  this  symposium,  for  some  further  comments 
on  the  need  for  an  adequately  fine  analysis. 
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distribution.  3ecause  in  practice  a  space  vehicle  is  likely  to  be 
drifting  or  slowly  tumbling,  the  results  computed  and  here  reported  are 
for  isotropic  distributions.  Also,  DOSE(S)  was  normalized  to  one 
proton/ca^  of  energy  S  >  100  Mev  (for  power  law  spectra),  so  the  normal¬ 
ization  of  F,  incident  flux,  must  be  compatibly  defined. 

Thus,  in  a  vehicle  shielding  design  problem,  a  geometrical  analysis 
is  made  of  the  proposed  design  to  determine  a  number  of  "constant¬ 
shielding"  regions  about  an  interior  point.  The  dose  at  that  point  is 
calculated,  one-dimensionally,  for  each  of  the  constant-shielding 
regions  and  the  weighted  sum,  Eq.  8,  is  performed  using  the  dose  values 
and  the  fractional  solid  angles  for  each  region.  The  vehicle  analysis 
part  is  a  drafting  board  problem,  and  can  be  done  quickly  and  with 
adequate  accuracy  for  most  practical  vehicle  configurations.  But  the 
power  of  this  method  is  in  the  fact  that  any  complex  arrangement  can  be 
analyzed  in  these  terms.  The  dose  results  presented  here  apply  to 

specific  points  inside  a  very  complex  shielding  arrangement — the  tissue 
3 

of  the  human  body. 

Secondary  Dose  Code 

Secondary  Radiation  Components.  In  addition  to  the  dose  due  to 
primary  penetrating  protons,  there  is  some  energy  deposition  by  the 
secondary  radiations  generated  by  the  incident  protons.  These  second¬ 
aries  include  gamma  rays,  and  evaporation  and  knock-on  nucleons.  The 
neutrons  and  protons  are  generated  in  nuclear  reactions  where  an  inci¬ 
dent  proton  excites  a  shield  material  nucleus.  Some  neutrons  and  pro¬ 
tons  of  high  energy  (up  to  about  the  incident  proton  energy)  con 
theoretically  be  produced.  The  bulk  of  the  secondary  neutrons  and 
protons,  perhaps  with  ten  times  the  intensity  of  the  high  energy  tail, 
are  produced  by  evaporation  from  the  excited  nucleus.  These  evaporation 
nucleons  are  generated  isotropically  and  have  energies  of  5  to  2$  Kev, 
and  proton  numbers  are  roughly  one  half  the  neutron  numbers.  The 
statistical  theory  of  the  nucleus,  which  formally  describes  the  exci¬ 
tation  and  evaporation  processes,  is  invalid  for  ths  low  atomic  number 
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nuclei,  ouch  as  carbon  or  aluminum,  so  the  proton  interaction  pro'oabili- 

7 

ties  were  not  computed  using  that  theory.  Rather,  experimental  data 
were  used  to  estimate  the  total  neutron  production  probabilities  as  a 
function  of  incident  proton  energy. 

Gamma  rays  emerge  from  excited  nuclei,  having  energies  depending 

on  the  specific  nuclei  and  excitation  levels.  The  excitation  may  be 

residual  from  a  higher  energy  nucleon-producing  interaction,  or  it  may 

be  from  an  inelastic  scattering  process  in  which  the  incident  proton 

loses  a  relatively  small  amount  of  energy  to  the  nucleus.  From  either 

cause,  a  discrete  gamma  ray  spectrum  results  from  the  passage  of  high 

8 

energy  protons  through  matter. 

The  secondary  dose  code  was  developed  to  compute,  for  all  three 
components,  the  secondary  fluxes  and  spectra  generated  in  successive 
layers  of  the  shielding,  the  secondary  flux  and  spectra  penetrating 
through  the  material  between  the  generation  layer  and  the  interior  point 
of  interest,  and  the  total  absorbed  dose  at  the  interior  point. 

Description  of  the  Code.  The  secondary  code  will  not  be  described 
here  in  detail,  as  it  is  to  be  reported  elsewhere.  Only  enough  dis¬ 
cussion  will  be  given  to  indicate  how  the  preliminary  results  depend 
on  the  assumed  secondary  radiation  production  data,  since  thfese  data 
are  likely,  to  be  better  known  as  more  nuclear  physics  research  is  done. 
Then,  presumably,  better  results  will  be  obtained. 


7.  W.  E.  Crandall  and  G.  P.  Millburn,  Neutron  Production  at  High 
Energies,  J.  Apol.  Phys.  29,698  (1958). 

See  also  R.  Wallace  and  C.  Sondhaua,  (especially  Figs.  12-16)  this 
symposium,  and  UCRL-10439* 

8.  R.  Madey,  A.  G.  Duneer,  Jr.,  and  T.  J.  Krieger,  Gamma  Dose  from 
Solar  Flare  Protons  Incident  on  an  Aluminum  Shield.  Trans.  Am.  Nuc. 
Soc.  June  1962  meeting. 

•  D.  L.  Dye,  Space  Proton  Doses  at  Points  Within  ths  Hunan  Body. 

The  Boeing  Company,  D2-90106  (1962). 
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S  DEGRADED  THROUGH  XS-X 
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Consider  a  layer  of  shielding  material,  Fig.  2,  of  total  thickness 
XS.  At  a  depth  X,  in  a  layer  dX  thick,  the  secondary  radiation  is  being 
generated  according  to: 


N(E  )  = 
n 


6  x  lO^dX 


/ P(Sp,X)  XSKUM(SpfSa)  Cra(Ep)  dEp, 


(9) 


23  /• 

r(Sr)  -  - ' jL~ — —  J  P(Ep,X)  XGMULT(Ey)  Of  (S  ) 'dE  ,  (10) 


and 


s  (e  )e 

p  s  s 


6  x  lO^dX 


I 


P(E  ,X)  XPMULT(E  ,E  )dS  CT  (E  )  dE  .  (11) 
p  p  s  s  a  p  p 


Where:  A  is  atomic  number  of  the  shield  material, 

N(E^)  is  the  flux  of  neutrons  generated  in  an  energy  group 

_ _  having  energy  Eq, 

P(Ej-)  is  the  flux  of  photons  generated  having  energy  E  , 

Sp(Ea)d£s  is  the  differential  energy  distribution  of  secondary 

protons  generated  having  energy  E  to  E  +  dE  , 

P(Ep,X)d Ep  is  the  proton  energy  spectrum  at  depth  X,  the  flux 

having  energy  between  E  and  E  +  dE  , 

P  P  P 

XNMULT(Ep,EQ)  ®’n(2p)  is  the  yield  of  neutrons  having  energy  Ea, 
per  incident  unit  flux  of  protons  having  energy  Ep, 
XGMULT(Ef)  Oy(E^)  is  the  yield  of  gamma  rays  having  energy  E^ 

(due  to  inelastic  scatter  or  other  processes)  produced 

XGM’JLT 

depends  primarily  on  the  nuclear  levels  of  the  struck 
nucleus,  while  Gy  is  the  cross  section  for  an  inter- 
action, 

XPMULT(E  ,E  )dE  O'  (E  )  is  the  differential  yield  of  secondary 

P  6  S  3  p 

protons  having  energy  E  to  E  +  dE  from  processes 
induced  by  incident  protons  of  energy  Ep. 

The  multiplicity  factors  (XNMULT,  XGMULT,  XPMULT)  are  normalized  to  one 
"particle"  produced  so  that  experimentally  determined  total  yield  cross 


per  incident  unit  proton  flux  of  Energy  Ep. 


sections  may  be  used. 


7,10 


10.  T.  Watasuki  et  al.,  Gamma  Rays  from  Several  Elements  Bombarded 
.  by  10  and  14  Mev  Protons,  J.  Phvs.  Soc.  Japan  1$,  1141  (I960). 
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A  geor.etri.oal  weighting  factor  (GECM)  is  assigned  to  the  radiation 
generated  in  any  given  material  layer.  This  factor  depends  upon  the 
angular  distributions  of  the  secondaries  and  upon  the  relative  sizes  of 
the  vehiole  shield  and  t'r.o  dose  point  of  interest.  The  fraction  of  the 
secondaries  produced  in  a  given  sublayer  within  the  shielding  that  is 
directed  toward  the  interior  dose  point  is  given  by  multiplying 
N  .  P  ,  or  S  by  GEQM  for  that  layer. 

O  5  p 

The  neutrons  and  gamma  rays  are  assumed  to  be  absorbed  by  inter¬ 
vening  shielding  according  to  exponential  laws.  The  protons  generated 

in  dX  are  added  to  P(E  ,X  +  dX)dE  .  (This  is  the  reason  we  needed  the 

P  P 

differential  energy  distribution  for  S  (E  )).  Thus,  the  spectrum  for 

•  P  P 

the  next  dX  is  enhanced  by  the  secondary  protons  generated  in  previous 
sublayers. 

In  the  actual  code  computation,  the  layer  XS  thick  is  broken  up 
into  NSL  finite-thickness  sublayers,  where  NSL  is  determined  by  a  speci¬ 
fied  maximum  sublayer  thickness,  except  that  NSL  never  exceeds  another 
specified  parameter,  NSUBL.  This  last  restriction  saves  computer  time 
for  thick  shields. 

The  absorption  of  the  neutrons  and  gammas  is  carried  through  all 
intervening  shielding  layers  to  the  dose  point  inside,  and  the  dose 
there  calculated  using  well-known  dose  conversion  factors.  The  code 
is  dimensioned  for  10  different  materials  and/or  layers.  The  secondary 
protons  are  also  carried  through  the  intervening  shielding  separately 
from  the  total  (primary  plus  secondary)  proton  flux,  so  as  to  show 
secondary  proton  dose  separately.  For  this,  the  proton  spectrum 
modification  and  dose  code  described  in  the  previous  section  is  used. 
After  the  absorbed  doses,  in  rsds,  are  computed  at  the  interior  point 
for  each  secondary  component,  these  doses  are  individually  converted 
to  biological  dose,  in  rem,  by  multiplying  by  the  appropriate,  specified, 
RBE.  The  individual  component,  as' well  as  the  total,  secondary  biodose 
is  then  read  out.  Finally,  the  weighted  sum,  Eq.  8,  is  done  to  obtain 
total  dose  from  all  directions  at  the  point. 
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Data  Used  lr,  the  Secondary  Dose  Computation.  In  addition  to  the 
proton  range-energy  tables  used  in  the  primary  dose  computation,  a 
number  of  nuclear  and  other  parameters  are  needed  in  the  secondary  code. 
These  include:  the  yield  cross  sections  O'  ,  0*^5  the  multiplicities 
XNKULT ,  XPMULT,  and  XGKULT;  the  discrete  gamma  ray  energies  generated 
in  the  specific  shielding  materials,  designated  E3AK;  the  gamma  and 
neutron  absorption  coefficients  and  dose  conversion  factors  for  these 
energies;  the  angular  distribution  factors,  GEOM,  for  each  layer;  and 
the  maximum  sublayer  thickness,  XLMAX ,  for  each  material.  Neutron 
yield  cross  sections  were  found  from  the  empirical  relation 


<r  a  C.E 
n  1 


,0.4 


n 

which  fits  the  experimental  data'  for  aluminum  and  carbon  (assumed  the 
same  as  tissue)  with  a  0.07  and  «*  0.035»  respectively.  (For 
iron,  a  0.3.)  Gamma  ray  yield  cross  sections  and  energies  were 
also  inferred  from  reported  experiments  and  the  Oy  values  used 

are  plotted,  along  with  0*n,  in  Fig.  3.  Two  different  fly  curves 
were  used,  one  in  which  cy  approaches  0*tota^  at  high  energies,  and 
one  in  which  it  is  held  constant  at  its  low  energy  maximum 
(  **  2  barns)  at  high  energies.  This  latter  case  is  labelled  "upper 
limit".  The  inelastic  scatter  gamma  ray  energies,  EG  AM,  and  multi¬ 
plicities,  XGKULT,  for  aluminum  and  tissue,  and  the  NBS  Handbook  dose 
conversion  coefficients,  are  given  in  Table  1. 


The  secondary  nucleon  energy  distributions  are  somewhat  less  well 
known,  and  as  better  data  become  available,  better  dose  values  can 
be  obtained  by  using  them  in  this  code.  For  the  present  calculation, 
an  evaporation  hump  ranging  from  5  to  23  Mev  for  secondaries  was 
assumed  for  all  incident  proton  energies  above  23  Kev,  and  a  high- 
energy  tail  was  assumed  to  extend  out  to  0.3  the  incident  proton  energy, 


11.  G.  Schrank  et  al.,  Inelastic  Scattering  of  17-Kev  Protons,  Phys. 
Rev.  127,  2139  (19^2).  This  paper  gives  many  references  to 
experiosntal  determinations. 
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Table  1.  Gamma  Say  Energies,  Multiplicities,  and  Dose  Coefficients 
Assumed  for  the  Secondary  Dose  Computation 


Material 

B3AK 

(Mev) 

XGMULT 

DOSGAM  (Tissue) 
(rad/proton/cm^ ) 

Aluminum 

1.4 

0.45  . 

4.8  x  10-9 

2.2 

0.19 

4.0  x  10“9 

3.0 

0.25  " 

3.5  x  10~9 

5.° 

0.08 

3-3  x  10"9 

r 

8.0 

0.05 

3.2  x  10-9 

Tissue 

4.4 

0.88 

3.0  x  10“9 

6.1 

0.06 

2.9  x  10‘9 

7.0 

0.06 

2.9  x  10“9 

coot 


The 


containing  about  10  percent  of  the  total  secondary  nucleons.* 
secondary  spectrum  was  normalized,  for  each  incident  proton  energy  in 
the  integration  of  Zq.  11,  to  one  nucleon  over  the  whole  energy  range 
1  ^  Ss  <  2^/2.  Mora  detail  is  given  elsewhere.^ 

The  GECM  factors  were  assumed  to  be  constant  for  each  shielding 
layer,  and  were  taken  as  0.5  for  both  outer  aluminum  and  inner  tissue 
layers  in  the  results  to  be  given  here. 

Hur.an  Geometry 

The  absorbed  dose  at  each  of  several  points  in  a  seated  astronaut  • 
is  desired.  The  points  of  specific  interest  are  those  that  are  parts 
of  radiosensitive  physiological  systems,  such  as  hematopoetic  and  spinal 
marrow,  gut,  and  eye  lens.  Twelve  such  body  points  were  chosen,  located 

as  indicated  in  Table  2,  on  a  75-percentile  seated  man.  This  body 

12 

geometry  is  based  on  data  given  in  an  anthropometric  study  of  Air  ?orce 
personnel.  (Seventy-five  percent  of  the  men  had  smaller  body  measure¬ 
ments  than  those  used.)  Scale  drawings  of  this  composite  astronaut's 
body  were  made  and  used  to  determine  the  tissue  thicknesses  around  each 
of  the  selected  body  points  in  all  directions."^  For  each  specific  body 
point  listed  in  Table  2,  this  descriptive  geometry  problem  was  solved; 
the  point  is  taken  as  the  center  of  a  large  sphere  on  which  isothickness 
contours  are  projected.  These  contours  mark  out  a  set  of  solid  angular 
regions;  a  region's  projected  area  on  the  sphere  divided  by  the  total 
sphere  area  is  the  fractional  solid  angle  of  that  region.  Depending 
on  the  accuracy  required,  one  may  choose  any  set  of  regions  to  define 
the  fractional  solid  angles  for  the  human  geometry  problem. 

.  The  isothickness  contours  are  obtained  by  drawing  a  series  of 
cross  sectional  views  of  the  body  at  different  angles  <f>  from  the 

12.  Hertzberg,  Daniels,  and  Churchill,  Anthropometry  of  Flying 
Personnel — 1950.  VADC  Technical  Report  52-321  (195^) • 

•The  form  of  this  secondary  nucleon  distribution  is  not  unlike  that 
given  in  this  symposium  by  Dr.  K.  Strauch  as  a  "hypothetical  particle 
spectrum",  built-up  on  the  high  energy  end  to  account  for  secondary 
angular  distributions. 


TABLE  2.  DESCRIPTION  OF  BODY  POINTS  ANALYZED 
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anterior-posterior  plane;  and  by  drawing  a  series  of  horizontal  cross 
sectional  views  at  various  planes  above  and  below  the  point  of  interest. 

A  tissue  thickness  is  then  obtained  by  determining  the  true  distance 
from  the  point  to  the  body  surface  along  every  (latitude)  direction  6 
in  the  plane.  Figure  4  illustrates  the  method  for  a  specific  point. 

The  result  of  this  descriptive  geometry  problem  is  a  map  on  the  pro¬ 
jection  sphere  showing  lines  of  constant  tissue  thickness  between  the 
point  of  interest  and  the  body  exterior  along  the  cf>  direction.  Figure  5 
shows  a  typical  isothickness  contour  for  a  body  point,  in  the  spinal 
column.  In  this  procedure,  the  human  body  is  assumed  to  be  a  homogeneous 
tissue  mass  except  for  certain  cases  below. 

Conceptually  simple,  but  computationally  tedious,  extensions  of 
the  interior  body  point  dose  calculations  for  the  STE2NUM,  FEMUR,  WAIST  4, 
and  EYE  LENS  positions  were  also  analyzed  with  slight  geometric  vari¬ 
ations.  These  variations  were  intended;  (1)  to  solve  more  realistic 
problems  that  include  bones  and  lung  absorption  effects,’  (2)  to  give 
an  idea  as  to  the  dose  error  made  in  neglecting  bones  and  lungs;  i.o., 
in  assuming  the  body  to  be  homogeneous,  (3)  to  assess  the  shielding 
effectiveness  of  partial  body  shielding  built  in  as  a  chair,  and  (4)  to 
assess  the  effect  of  goggles  in  eye  protection. 

\ 

Incident  Proton  Spectra  » 

Five  different  incident  proton  spectra  were  considered;  four  having 

power  law  forms,  the  other  approximating  the  Van  Allen  belt  spectrum. 

1  o 

The  power  law  forms,  normalized  to  one  protoV®®  of  E  >  100  Mev,  are 

given  by: 

N(E)dLE  -  <n  -  l)(100)a  "  1  E“n  d£  for  20  Mev  <,  E  £  1,000  Kav 

for  the  parameter  n  with  values  2,  3,  4,  and  3*  They  are  flat  at  the 

20  Mev  value  for  the  interval  1  ^  E  £  20  Mev.  The  Van  Allen  spectrum 

2  *2 

is  normalized  to  one  proton/cm  of  E  >  40  Mev,  and  has  an  E  form 
above  100  Mev,  with  the  lower  energy  part  of  Hess's  theoretical  form 
matched  to  it  at  100  Mev. 


■>  FRONT 


SIDE  ViEW,  p  =  90e 


"GUT"/' 
a.  POINT 


HORIZONTAL 

PLANE 


£  IS  VERTICAL  AXIS,  PASSING  THROUGH  THE  POINT  OF  INTEREST 

p  IS  LONGITUDE  ANGLE  - 

©  IS  LATITUDE  ANGLE 

X  ‘IS  THICKNESS,  A  FUNCTION  OF  0  AND  0  FOR  THIS  POINT 


FIGURE  4 
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pQ33  Results 

Figure  6  shows  the  penetrating  primary  proton  dose  and  the  secondary 
dose  plotted  separately  for  the  E-'*  spectrum  and  for  one-dimensional 
cases.  It  is  seen  that  the  nuclear  data  assumptions  made  here  lead  to 
a  relatively  small  secondary  radiation  dose,  even  with  "upper  limit" 
inelastic  scatter  cross  sections.  It  should  be  noted  that  for  steeper 
primary  spectra  the  primary  absorption  is  greater,  so  that  the  secondary 

dose  is  a  larger  fraction  of  the  primary  than  for  spectra  having  smaller 

9 

n  values.  However,  the  secondary  dose  is  clearly  a  small  fraction  of 

2 

the  primary  proton  dose  for  tissue  thicknesses  up  to  20  gm/cm  and 
aluminum  shields  of  up  to  10  gq/cm^,  for  typical  space  proton  spectra. 
Further,  the  solid  angle  subtended  by  tissue  thickness  greater  than 
20  gm/cm^  around  the  deepest  body  points  considered  (GUT,  WAIST  8)  totals 
only  about  20  percent^,  and  it  is  even  less  for  most  radiosensitive 
organs.  Thus  the  total  dose  to  these  organs  is  essentially  equal  to  the 
primary  proton  dose,  computed  by  Eq.  8,  for  feasibly  thin  external 
aluminum  shields. 

Tables  3  through  7  show  the  total  doses  at  each  of  the  body  points 

-2  -3  -4  -5 

considered  for  a  variety  of  aluminum  shells  for  the  3  ,  E  ,  E  ,  E  , 

and  Van  Allen  spectra,  respectively.  The  units  are  tissue  rad  per 
incident  isotropic  proton/cm^  having  E  >  100  Kev.  These  valuqs  utilize 
the  three-dimensional  human  geometry  analyses,  Eq.  8,  and  the  solid¬ 
line  cross  section  data  of  Fig.  3.  The  nuclear  data  used  are  thus 
reasonable  composites  of  published  experimental  information. 

In  addition  to  the  twelve  body  points  with  only  external  shell 
shielding,  four  special  cases  are  listed.  In  the  case  WAIST  4S,  a 
shielding  chair  of  either  2  gq/cm^  or  4  gq/sm^  polyethylene  is  con¬ 
sidered  with  three  outer  aluminum  shell  values.  The  line  labelled 
FS34UR  B  is  a  72-region  case  in  which  the  structure  of  the  bone  was 
considered.  The  line  labelled  EYE  G  is  a  case  in  which  0.23  cm  of 
lead  glass  was  added  to  the  zero  tissue  thickness  region  of  the  EYE 
analysis;  that  is,  to  the  front  of  the  eye  lens.  The  case  CHEST  2 
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PRIMARY  PROTON  DOSE  vs  LATERAL  WAIST  DEPTH  FOR 
VARIOUS  ALUMINUM  SHELL  SHIELD  THICKNESS.  E3  SPECTRUM 
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ignoren  air  in  the  lungs  and  bilateral  asymmetry  due  to  the  heart; 
CHEST  2L  takes  these  into  account.  It  is  aeon  fro.?.  Tables  3  through  7 
that  to  take  account  of  lung  air  or  bone  structure  is  hardly  worth  the 
extra  analysis  effort  and  computing  time  for  the  proton  dose. 

The  depth-dose  pattern  in  the  human  body  is  felt  to  be  of  special 
significance"'"'*  and  so  the  doses  at  points  laterally  along  the  waist 
are  plotted  in  Fig.  ?  to  show  the  representative  pattern.  It  is  seen 
from  Fig.  7  and  the  tabulated  values  that  a  midline  dose  of  15  to  30 
percent  the  skin  dose  is  typical,  dpendir.g  on  shielding  and  spectrum 
in  obvious  ways. 

Because  of  the  normalization  used  in  the  values  reported  here, 

they  may  readily  be  applied  to  environment  data.  If  a  trajectory  is 

known,  the  trapped  proton  flux  of  E  >  40  Kev  may  be  computed  and  inte- 

14 

grated  over  the  orbit.  If  a  specific  solar  flare  proton"  event  is 
estimated  to  have  a  spectrum  denoted  by  n  Bind  size  in  protons/cm^, 
both  parameters  integrated  over  the  whole  event,  (or  instantaneous 
flux  in  protons/cm  -sec),  then  total  dose  (or  dose  rate)  may  be  deter¬ 
mined  to  any  of  the  oody  points.  An  estimated  flare  size-probability 
curve  can  be  converted  to  significant  dose-probability^  using  these 
data  and  vehicle  shielding  assumptions.  By  simple  descriptive  geometry 
techniques  and  the  codes  here  reported,  the  shielding  required  in 
specific  vehicle  designs  may  be  computed  so  as  to  maintain  a  desired 
lower  probability  limit  for  specified"  organ  tolerance  dose. 
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